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Abstract. We propose two viable scenarios explaining the recent observations on cosmic 
positron excess. In both scenarios, the present relic density in the Universe is assumed to 
be still supported by thermally produced WIMP or LSP (x)- One of the scenarios is based on 
two dark matter (DM) components (Xi^O scenario, and the other is on SO(10) SUSY GUT. 
In the two DM components scenario, extremely small amount of non-thermally produced meta- 
stable DM component [0(1O~ 10 ) < nx /n x ] explains the cosmic positron excess. In the SO(10) 
model, extremely small R-parity violation for LSP decay to e is naturally achieved with a 
non-zero VEV of v c and a global symmetry. 



1. Introduction 

For a long time thermally produced weakly interacting massive particles (WIMPs) have been 
believed to be the most promising dark matter (DM) candidates. It is because the correct order 
of the magnitude of the cross section for explaining the present relic density of the Universe is 
naturally possible with WIMPs. Particularly, the lightest supersymmetric particle (LSP), which 
is a well-motivated particle coming from the promising particle physics model, i.e. the minimal 
supersymmetric standard model (MSSM), has attracted much attentions as an excellent example 
of WIMP. Actually, the Universe relic density by WIMP or LSP DM is the traditional scenario, 
which has been believed so far. 

However, recently some astrophysical experimental groups including PAMELA [3], ATIC 
[1], and the Fermi-LAT collaborations [5] reported the very challenging observations in cosmic 
ray: PAMELA observed positron fractions [e + /(e + + e - )] exceeding the theoretical expectation 
[6] above 10 GeV upto 100 GeV. On the other hand, the PAMELA'S observations on anti- 
proton/proton flux ratio were quite consistent with the theoretical calculation. The ATIC and 
Fermi-LAT's observations exhibit excesses of (e + + e~) flux in cosmic ray from 100 GeV to 1 
TeV. They would result from the positron flux that keeps rising upto 1 TeV. 

Apparently the above observational results are very hard to be interpreted in view of the 
conventional MSSM cold DM scenario: explaining the excess positrons with annihilations of 
Majorana fermions such as the LSP needs a too huge boost factor. Moreover, ATIC and Fermi- 
LAT's observations seem to require a TeV scale DM, if they are caused indeed by DM annihilation 
or decay. However, TeV scale DM seems to be disfavored by the gamma ray data, if the excess 
positron flux should be originated from the same physics explaining DM creation in the early 
Universe, i.e. from DM annihilation [TJ. 



1 Talks presented in PASCOS, SUSY, and COSMO/CosPA in 2010, based on Refs. Q][2]. 



Such astrophysical observations seem to destroy our traditional scenario on DM. In this 
article, however, it will be pointed out that such observations do not necessarily imply the 
presence of a new DM theory replacing the traditional one. We will propose two viable 
scenarios explaining them, in which the present relic density in the Universe is assumed to 
be still supported by thermally produced WIMP. 



2. Two dark matter components scenario 

Let us consider DM decay scenario to explain the cosmic positron excess. "Helicity suppression" 
is not valid in DM decay any longer. Unlike the DM annihilation scenario, DM decay scenario is 
relatively free from the gamma ray constraint, since the positron flux is just linearly proportional 
to the number density of DM [8]. In the DM decay scenario, however, there are some serious 
hurdles to overcome: (1) one is to naturally obtain the extremely small decay rate of the DM, 
Tdm ~ 10 -26 sec. -1 , and (2) the other is to naturally explain the relic density of the DM. 

The first hurdle could be somehow resolved by introducing an extra symmetry, an extra DM 
component with a TeV scale mass, and grand unified theory (GUT) scale superheavy particles, 
which mediate DM decay into the standard model (SM) charged leptons (and the LSP) [2]. It 
is because the required decay rate of 10 -26 sec. -1 can be achieved, if the dominant operator for 
decay to e + e - is dimension 6 suppressed by Mq UT , i.e. a four fermion interaction: 

5 

T m DM ln -26 ^ -1 m 

FdM ~ 196vr3il4 UT ~ 10 SeC ' ' (1) 

where ttidm ~ a few TeV. The fact that the GUT scale particles are involved in the DM decay 
might be an important hint supporting GUT [9]. However, since the interaction between the new 
DM and the SM charged lepton are made extremely weak by introducing superheavy particles 
mediating the DM decay, it is hard to thermally produce the new DM, and so non-thermal 
production of DM with a carefully tuned reheating temperature should be necessarily assumed 
for the required relic density. One way to avoid it is to consider a SUSY model with two DM 
components (%, X), where x 1S j us t the ordinary WIMP such as the LSP and X indicates a new 
DM component [2 [9]: 

X '• main component of DM explaining the relic density, thermally produced, 

absolutely stable. 
X : (extremely) small number density, non-thermally produced, 

meta-stable, decay to e + e explaining PAMELA/Fermi-LAT. 

In this class of models, a (global) symmetry should be introduced to forbid all unwanted 
dimension 4 and 5 operators contributing to DM decay. The required number density of X 
(= nx) turns out to be just 0(1CU 10 ) < nx/n x . It is possible only if the suppression of the 
dimension 6 operator is smaller than Mgut> 10 12 GeV < < 10 16 GeV [2], since the 

positron flux in the case of DM decay is proportional to nx ■ ^x- The low energy field spectrum 
in this class of models [2j is the same as that of the MSSM except for the neutral singlet extra 
DM component. Moreover, the models of [2] can be embedded in the flipped SU(5) GUT. 



3. LSP DM scenario: SO(10) SUSY GUT model 

In the second scenario, we suppose again that just the conventional bino-like LSP is the main 
component of the DM. Since the "bino" is a WIMP, thermally produced binos could explain well 
the relic density of the Universe. Without introducing a new DM component and interaction, we 
will attempt to explain the PAMELA'S observation within the framework of the already existing 
particle physics model, SO(10) SUSY GUT. 



3.1. SU(5) vs. SU(2) R scales 

In terms of the SM's quantum numbers, the SO(10) generator (= 45^) is split into the 
generators of the SM gauge group plus {{l,l)_ x ,(l,l) x }, (1,1) , {(3, 2)_ 5/6 , (3, 2) 5/6 }, and 
{(3, 2) X / G , (3, 2)__ 1 y 6 ; (3, l) 2 /3> (3, l)_ 2 / 3 }. We will simply write them as 

{E,E C }, N, {Q',Q' C }, {Q,Q c ; U,U C } , (2) 

respectively. The SM gauge group's generators and {E,E C }, N compose the generators of 
SU(3) c xSU(2) L xSU(2) fi xU(l)B_ L (= LR), where {E,E C } and a linear combination of the 
SM hypercharge generator and N (= Nr) correspond to the SU(2)/j generators. The other 
combination orthogonal to it is the \J(1)b—l generator (= Nbl)- Note that {E,E C } and N 
don't carry color charges. By the VEV of the adjoint Higgs (45//), the SO(10) gauge symmetry 
may break to LR. On the other hand, the SM gauge group's generators and {Q', Q c } compose 
the SU(5) generators. The VEVs of the Higgs in the spinorial representations (16//), (16//) 
breaks SO(10) down to SU(5). Hence, (45//) and (16//) determine the SU(5) and LR breaking 
scales, respectively. 

Non-zero VEVs of both (45//) and (16//) eventually give the SM gauge symmetry at low 
energies! Thus, if (45//) > (16//) = (16//) / ((16//) = (16//) > (45//) / 0), the gauge 
bosons and gauginos of {Q' , Q c } achieve heavier (lighter) masses than those of {E, E c } and N. 
The masses of the gauge sectors for {Q, Q c ; U, U c } would be given dominantly by the heavier 
masses in any cases, since both (45//) and {(16//), (16//)} contribute to their masses. 

3.2. LSP decay in SO (10) 

To achieve the needed extremely small decay rate of the bino-like LSP Xi we need extremely 
small R-parity violation naturally, x can never decay, if (1) R-parity is absolutely preserved and 
(2) x is really the LSP. We mildly relax these two conditions such that x can decay by assuming 

(1) a non-zero VEV of the superpartner of one family of RH neutrino, v x 
(i.e. R-parity violation), or 

(2) a mass of v x lighter than the x' s mass, m x (i.e. v x LSP). 

By introducing a global symmetry, one can forbid the (renormalizable) Yukawa couplings 
between v x and the MSSM fields. Then, v x can interact with the MS SM fields only through the 
superheavy gauge fields and gauginos of SO(10), since the (s)RH neutrino u x (p x ) is a neutral 
singlet under the SM gauge symmetry but it is charged under SO(10). It is embedded e.g. in 16 
of SO(10). Consequently, decay of x is possible but extremely suppressed. For instance, refer to 
the diagram of Figure [!}■ (a). We will discuss how this diagram can be dominant for the x decay. 

3. 3. The conditions for leptonic decay of x 

Let us consider the interactions of the superheavy gauginosfl In Table Q3 we list all the gauge 
interactions between the superheavy gauginos of SO(10) and two MSSM fields. They are, of 
course, the renormalizable operators. As seen in Table Q3 uf or uf couples to the superheavy 
SO(10) gauginos, {E,E C }, N, {Q,Q C }, and {U,U C }. 

According to PAMELA data [3] , the branching ratio of the hadronic DM decay modes should 
not exceed 10 %. To make the leptonic interactions, i.e. e^*vfE c , i/f*efE, and vf*vfN, e c *e\N 
dominant over the other interactions in Table [TJ we assume that 

2 Alternatively, one can employ the large representations, 126h, 126h, and 210h, instead of 16h, 16.tr, and 
45 H . 126 h and 126 H break SO(10) to SU(5), while 210 H breaks SO(10) to SU(4) C xSU(2) L xSU(2) ii . In our 
discussion throughout this article, 16jj (16 h) and 45h can be replaced by 126^ (126h) and 210.fr, respectively. 

3 In this article we don't discuss the cases in which x decays through the mediation of the superheavy gauge 
bosons. However, it turns out that the decay channels of \ through the mediation of the superheavy gauge fields 
are relatively suppressed, compared to the mediation of the superheavy gauginos discussed here [I]. 



(a) (b) 

Figure 1. Dominant diagram of the bino decay (a) and the gauge interaction between 
electrically charged superheavy LR gauginos and the MSSM lepton singlets (b). 
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Table 1. Gauge interactions between two MSSM fields and a heavy gaugino in the SO(10) GUT 



• The LR (or B — L) breaking scale should be lower than the SU(5) breaking scale, i.e. 
(16//) "C (45//). Then Mqi, Mq, Mrjj| become much heavier than Me and M/v, and so most of 
hadronic decay modes of x can be easily suppressed except those by E c , E, and N in Table [TJ 

• The slepton ef , which composes an SU(2)# doublet together with uf, needs to be lighter 
than the squarks. Then the decay channels of x by df*ufE c , uf*dfE, and uf*ufN ', d c *d\N become 
suppressed. We also require that x an d ef are much lighter than the charged MSSM Higgs. So 
the leptonic interactions, ef*ufE c , z/fef-E 1 , and uf*vfN, efefiV can dominate over the others. 

• At least one RH neutrino, i.e. the SU(2) l singlet neutrino uf (and its superpartner uf) must 
be lighter than x so that X decays to charged leptons. It is because uf is always accompanied 
by uf in the effective operators leading to the leptonic decay of Xi composed of efufE , uf*e\E, 
and uf*ufN, ef* ef N. If all the sneutrino masses are heavier than x-, v \ must develop a VEV for 
decay of x- Once uf is light enough, uf can achieve a VEV much easily. 

To be consistent with PAMELA'S observations on high energy galactic positron excess [3], 
the DM mass should be around 300 - 400 GeV [TT]! Thus, we simply take the following values; 

4 The triplets Higgs contained in 10^ can achieve the mass proportional to (45.tr) via 10^45^10^, |10j . 

5 The (e + + e~) excess observed by Fermi-LAT could be explained by astrophysical sources such as nearby 
pulsars (and/or with the sub-dominant extra TeV scale DM component [2]). In fact, pulsars can explain both the 
PAMELA and Fermi-LAT's data in a suitable parameter range 12 . However, this does not imply that DM in 



1. (16//) < (45//). (If m x < mpc, then (i/f) ^ 0.) We will assume (16//) ~ 10 14 GeV. 

2. squarks, charged Higgs, higgsinos and other typical soft masses are of 0(1) TeV. 

3. m u c < m x ~ 300 - 400 GeV < mgj « 0(1) TeV. 

Consequently, SO (10) is broken first to LR, which would be the effective gauge symmetry valid 
below the GUT scale. As seen from Table [TJ the gauge interactions by the LR gauginos (and 
also gauge fields) preserve the baryon numbers. Even if the masses of the LR gauginos and 
gauge fields are relatively light, their gauge interactions don't give rise to proton decay. 



3-4- Seesaw mechanism 

Although one RH neutrino is light enough, the seesaw mechanism for obtaining the three 
extremely light physical neutrinos still may work. Let us consider the following superpotential; 

W v = y$ kK^(j / 1) + l -M l3 ufu^i, j^l), (3) 

where the Majorana mass term of i/f could be generated from the non-renormalizable 
superpotential (16#)(16#)16jl6j/Mp ^ 1). Thus, My (3> {h u }) could be determined, 
if the LR breaking scale by (16//) is known. In this superpotential, we note that one RH 
neutrino v\ does not couple to the MSSM lepton doublets and Higgs. For instance, by assigning 
an exotic U(l) R-charge to z/f, one can forbid its Yukawa couplings to the MSSM superfields. 
Thus, v\ would be decoupled from the other MSSM fields, were it not for the heavy gauge fields 
and gauginos of the SO(10) SUSY GUT. 

Taking into account only Eq. @, one neutrino remains massless. The two heavy Majorana 
mass terms of and are sufficient for the other two neutrinos to achieve extremely small 
physical masses through the constrained seesaw mechanism [13J: 



m v = mE = - I v 2 2 v 2 3 



Vn Vl 3 
V22 V23 
t>32 V33 



/ \ / 

M£ M^ 1 v 12 v 2 2 v 32 ] , (4) 

M^ 1 M33 1 / V "13 V 23 V 33 



where Vij = (h u ), and M- 1 denotes the inverse matrix of Mij. One of the eigenvalues of m v 
is zero and the other two are of order v 2 /M. Through the diagonalization of the mass matrix in 
Eq. (jlj), the three left-handed neutrinos from the lepton doublet l±, I2, and I3 can be maximally 
mixed, whereas the mixing of the RH neutrinos is only between and z/|. A complex phase in 

yf^ could make the leptogenesis possible [13]. 



3.5. LSP decay rate and the seesaw scale 

Let us consider the following terms in the superpotential; 

-^-(Te^ie^ + KE 3 , (5) 

where Mp = 2.4 x 10 18 GeV and k is a dimensionless coupling constant. S is an SO(10) singlet. 
We assign e.g. the U(l) R-charge of 2/3 to I61 and S, and to 16//. This charge assignment 

addition to pulsars can not be the source of the galactic positrons. In fact, we don't know yet a complete pulsar 
model, in which all the free parameters would be fixed by the fundamental physical constants. 

Alternatively, one could assume m x ~ 3.5 TeV and the model is slightly modified such that \ decays dominantly 
to fj, , 1/2 rather than to e , uf, which is straightforward, the Fermi-LAT's data as well as the PAMELA'S can 
be also explained [TT]. In this case, however, the soft SUSY breaking scale should be higher than 3.5 TeV. 



forbids the renormalizable Yukawa couplings between v\ and other MSSM fields carrying integer 
R-charges. 

The scale of (16//) (= (16//) = Me/ V2gio) can be determined such that it is consistent with 
PAMELA data. The soft mass term of E and the A-term of kS 3 in the scalar potential permit 
a VEV (£) ~ "1,3/2/ K - Then, the scalar potential generates a linear term of i>\ coming from the 
A-term corresponding to the first term of Eq. ([5]), V D mjy 2 ((16//)/K 2 Mp)z}]\ The linear term 
and the soft mass term of v\ in the scalar potential can induce a non-zero VEV of v\: 



^ ^ X Mp- (6) 



Thus, the decay rate of x i n Figure CD- (a) can be estimated: 
r _ a\ a Y m x I m 3/2 (^) \ 2 _ a 2 w a Y m x ( m 3/2 



~ 10~ 26 sec." 1 , (7) 



x 96M| I m% I mM E M 2 p 



Kmzc 
1 , 



where aio (= g 2 /47r) and ay [= g Y /^ = (3/5) x gl/Air, where g\ is the SO(10) normalized 
gauge coupling of g Y ] are approximately 1/24 and 1/100, respectively. Here, we ignore the 
RG correction to a\$. 300 - 400 GeV fermionic DM decaying to and a light neutral 
particle can fit the PAMELA data [llj. For m x rj 300 - 400 GeV, {m^^j Km- e () ~ 10, 
Me or (16//) is estimated to be of order 10 14 GeV. This is consistent with the assumption 
(16//) < (45//) ~ 10 16 GeV. Therefore, the masses of the other two heavy RH neutrinos, which 
do not contribute to the process of Figure [1]- (a), are around 10 10 GeV or smaller in this case: 
W D ^•((16//)(16//)/M P )16 i 16 i (^i ^ 1) D ^(10 10 GeV) x v^(i,j + 1). 

4. Conclusion 

In this article, we pointed out that the traditional DM scenario based on thermally produced 
WIMP do not necessarily conflict with the cosmic positron excess observed by PAMELA and 
Fermi-LAT. We have proposed two viable scenarios based on two DM components and SO(10) 
SUSY GUT. Particularly in the SO(10) model, extremely small R-parity violation for LSP decay 
could be naturally achieved with a non-zero VEV of v\ and a global symmetry. 
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